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1 INTRODUCTION 

The gas of the interacted electrons is usually described within Kohn-Sham approximation by 
the set of Poisson and Schrodinger equations with an effective potential for the single-particle 
wave functions. The solution can be obtained using many-step iteration procedure. The well 
known difficulty in this task is that the wave functions obtained after every iteration step give 
the distribution of electron density which is not correspond to the boundary conditions for 
the Coulomb potential. As a result, either it is impossible to obtain the solution for the next 
iteration step or some parameters of the system are to be changed, for example, the density 
of the positive charge. The proposed new iterative scheme for solving Kohn-Sham equations is 
freed up of necessity to modify parameters of the system while iteration process. 

2 SELF-CONSISTENCY PROCEDURE 

The self-consistent electron density is constructed as the solution of the Kohn-Sham equations. 

^VVe(^) + iE- u,siz))Mz) = (1) 

\/\ = (Nn - niz)) (2) 

Here A^d is the density of the positive background, u is the Coulomb potential energy of electron 
and Uxc is the exchange-correlation potential energy, which we assume in the local-density 
approximation. 

We represent the total electron density as a sum 

n{z) = Uindiz) + n^niz), (3) 
where induced electron density depends on the potential as follows 
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This expression is the usual quasi-classical solution of Eq.(|I|) averaged over Fermi wavelength 
scale and 

E^Ef 

will be named as the quantum electron density. 

Using Eq.(^ Poisson equation can be rewritten in a new form 

V\ + nindM = A/'d - n^uiz) (4) 
The total iteration scheme can be described as following 

^ = 0,1,..., < = 

IV'^Eiz) + (E - ul^izWEiz) = 
^qu^ = Ee^Ep li'hWz) - n[^^ [uIq] 

This procedure was tested in solving the self-consistent set of the equations for two tasks within 
the jellium model: 

1. Semi-infinite electron gas which is bounded by self-consistent potential 

2. Semi-infinite electron gas which is bounded by infinite potential barrier. 



3 SEMI-INFINITE ELECTRON GAS BOUNDED BY 
SELF-CONSISTENT POTENTIAL 

Using the local-density approximation for exchange and correlation potentials we calculated 
the self-consistent distribution of the electrons and the work function in the wide range of the 
bulk electron densities. All computations were done in the framework of the jellium model, 
where the positive ions are replaced by the uniform background of positive-charge density. 
This background fills the half-space z > and is neutralized by the electron density. Figure 
|l] shows the distributions of electron density and effective potential calculated self-consistently 
for Rs = 4. These curves are compared with the corresponding results from where the 
straight iteration scheme was claimed not converging (see |jl|], Appendix B) and authors did 
have to use, in fact, variational method for solving Eqs. (^^. Our iteration scheme @ allows 
to obtain the self-consistent answer through the straight solution of Kohn-Sham equation set. 
According to the well-known Budd-Vannimenus theorem [@], the self-consistent potential must 
satisfy to the follow sum rule 

Ab = m(0) - u(oo) = n^^. (6) 

dn 

Here e is an energy per one electron for the homogeneous electron gas. Substituting the explicit 
expression for e, it is easy to calculate Ab as the function of the electron density. The same 
parameter obtained from our self-consistent solution we denote as Age- As it can be seen in 
Table I, parameters Ab and Age are in a good agreement. This fact is the independent evidence 
that our solution is fully self-consistent |Q| . The work function is defined as 



W = Meflf(-Oo) - /X. 
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Comparative plot of distributions of electron density and effective potential 

near the surface of simple metal, r^=4 
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Figure 1: 



The dependence of the self-consistently calculated Wsc on Wigner-Seitz parameter Rg is shown 
in Figure H in companion with results of not fully self-consistent calculations in This 
figure provides comparative analysis of work function calculations obtained by different tech- 
nique. 



Compartive plot of density dependence 
of work function within jellium model 
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TABLE I. Electron density dependence of the work function 





Wsc(eV) 


Ab 


Asc 


0.5 


3.21 


0.358 


0.358 


1.0 


3.40 


0.317 


0.317 


1.3 


3.51 


0.291 


0.291 


1.5 


3.60 


0.274 


0.274 


1.65 


3.63 


0.261 


0.261 


1.8 


3.64 


0.248 


0.248 


2.07 


3.60 


0.224 


0.224 


2.3 


3.53 


0.204 


0.204 


3.28 


3.12 


0.115 


0.115 


3.99 


2.87 


0.051 


0.063 


4.96 


2.51 


-0.048 


-0.041 


5.23 


2.46 


-0.075 


-0.074 



4 SEMI-INFINITE ELECTRON GAS BOUNDED BY 
INFINITE POTENTIAL BARRIER 

In this section we apply our iteration scheme for calculating self-consistent electron density 
in the case of electron gas bounded by infinite potential barrier. The positive charge and the 
electrons fill the half-space z > 0. The task was considered with accounting for the external 
electric field perpendicular to the surface plane. The presence of the electric field modifies the 



boundary condition for Coulomb potential at the surface. Using the results of this task, the 
capacitance of a tunnel structure can be estimated. Figure ^ demonstrates the distributions of 



^ Distributions of electron density, Coulomb and effective potentials 

^ near infinitely high barrier in an external electric field, R =2.07 
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Figure 3: 

electron density, effective and Coulomb potentials for Rg = 2.07 calculated by our technique. 
The solution of this task can be used for an estimation of the capacitance of a real finite-height 
barrier structure. The structure consists of two n— doped semiconductor half-spaces, separated 
by a potential barrier of width d. The total capacitance is C = dQ/dV, where Q is a total charge 
in the semiconductor lead and V = u{oo) — m(— cxd) is the voltage drop over the structure. Our 
calculations corresponds to the non-transparent infinitely high barrier with d = and give us 
the maximum achievable specific capacitance in barrier structures. In Figure ^ we plot the 
dependences of the capacitance on i^^calculated within the different approximations. 
The partial support from the Russian Foundation for Basic Researches is gratefully acknowl- 
edged. 



Density dependences of the capacitance calculated within 
different approximations 
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Figure 4: 
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